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BASED RANDOM NUMBER GENERATORS 

RELATED APPLICATIONS 
5 The following applications of the common assignee, which are hereby 

incorporated by reference, may contain some common disclosure and may relate to 
the present invention: 

U.S. Patent Application Serial No. / , entitled "RANDOM 

NUMBER GENERATORS IMPLEMENTED WITH CELLULAR AUTOMATA" 
fjj 10 (Attorney Docket No. 10017475-1); and 

U.S. Patent Application Serial No. / , entitled "SOFTWARE 

IMPLEMENTATION OF CELLULAR AUTOMATA BASED RANDOM 
NUMBER GENERATORS" (Attorney Docket No. 10019023-1). 
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1 5 FIELD OF THE INVENTION 
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This invention relates generally to random number generation. More 
specifically, this invention relates to systems and methods to generate cellular 
automata based random number generators (CA-based RNG). 

20 BACKGROUND OF THE INVENTION 

Since the inception of computers, random numbers have played important 
roles in areas such as Monte Carlo simulations, probabilistic computing methods 
(simulated annealing, genetic algorithms, neural networks, and the like), computer- 
based gaming, and very large scale integration (VLSI) chip-testing. The bulk of the 

25 investigation into random (more properly, pseudo-random) number generation 
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methods has been centered around arithmetic algorithms. This is because the 
prevalent computing medium has been the general purpose, arithmetic computer. 
Digital hardware designers have long relied on feedback shift registers to generate 
random numbers. 

With the advent of VLSI design, built-in self-tests have become advantageous. 
In this design, the bulk of the chip testing system is incorporated on the chip itself. 
Linear feedback shift registers were used initially to implement the random number 
generation portion of the built-in self-test. 

In 1986, Wolfram (S. Wolfram, "Random sequence generation by cellular 
automata," Advances in Applied Mathematics, vol. 7, pp. 123-169, June 1986) 
described a random sequence generation by a simple one-dimensional (1-d) cellular 
automata with a neighborhood size of three. The work focused on the properties of a 
particular CA-based RNG dubbed "CA30," so named due to the decimal value of its 
truth table. Statistical tests indicated that the CA30 was a superior random number 
generator to the ones based on linear feedback shift registers. Wolfram suggested that 
efficient hardware implementation of the CA30 should be possible. 

Hortensius et al. (P. D. Hortensius, R. D. McLeod, and H. C. Card, "Parallel 
number generation for VLSI systems using cellular automata," IEEE Transactions on 
Computers, vol. 38, no. 10, pp. 1466-1473, October 1989) described the use of CA30 
as a random number generator in an Ising computer. They also described using 
combinations of CAs (CA90 and CA150), which generated even better random 
numbers than the CA30. They further indicated that time and site spacing may 
improve statistical quality of random numbers generated by the CA. Time spacing is 
where the RNG is advanced more than one step between random number samples and 
site spacing is where not every bit value generated is used. 

i 
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Cellular automata (CA) may be thought of as a dynamic system discrete in 
both time and space. CA may be implemented as an array of cells with homogeneous 
functionality constrained to a regular lattice of some dimension. For example, in one- 
dimension, the lattice could be a string (open-ended) or a ring (close-ended), or in 
two-dimensions, the lattice could be a plane (open-ended) or a toroid (close-ended). 
Open-ended CAs have boundaries that are Fixed and close-ended CAs have 
boundaries that are periodic. 

A function of a CA cell may be represented as a truth table. Figure 1 A shows 
an exemplary truth table for a four-input CA cell. Figure IB shows an exemplary 
implementation of a cell of the CA. As shown, the cell i implicitly includes a one-bit 
register. In this instance, there are 1 6 possible conditions to which a cell may respond 
(the neighborhood size N is 4 corresponding to the number of inputs). The number of 

unique responses is 2^ N or 65,536 (see Table 1 below). In other words, there can be 

65,536 unique four-input machines for a given interconnection topology. 

Referring back to Figure 1A, a notation is provided to identify the CA 
implementing the above function. In essence, the output of the truth table is used as 
the identification in conjunction with the interconnection notation. As shown, the 
output of the truth table is converted to a number (from binary to base 16 to decimal). 
The CA represented by the truth table in Figure 1 A is denoted to be CA06990. 

As indicated before, a CA may be made of multiple cells, and the inputs of 
one cell may connected to the output of other cells. There may even be a feedback 
contact meaning that one of the inputs of the cell is connected the output of the cell 
itself. Thus, to uniquely identify a CA, the interconnection topology information 
should also be provided in addition to it's truth table representation. Figure 1C 
illustrates an exemplary notation, a relative displacement notation, which indicates the 
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interconnection topology information of cell i, i.e., how far away the connecting cells 
are relative to a given four-input cell i. 

As an example, Figure ID illustrates a 64-cell one-dimensional ring automata 
network with a displacement of {-1, 0, 1, 2} from the perspective of cell 0. In this 
instance, each cell i is assumed to have the same displacement value, i.e., all cells 
have identical functions. In a one-dimensional ring CA network, each cell i has two 
adjacent neighbors, one on either side. Because the CA network is periodic, cell 63 is 
adjacent to the cell 0, and thus the displacement of i - 1 from cell 0 lands on cell 63. 

In a one-dimensional CA network, a relative displacement value {-1, 0, 1, 2} 
indicates that dg input of cell i is connected to the output of the cell i - 1 (one cell to 
the left), the 04 input is connected to the output of the cell i itself, the d 2 input to cell 
i + 1, and the di input to cell i + 2. More specifically, from the perspective of cell 0, 
the inputs d 8 , 04, d 2 , and di are connected to the outputs of cell 63, itself, cell 1, and 
cell 2, respectively. 

Each cell in the CA network has a state that is updated as a function of its 
neighbor connections at each time step. In other words, the state of a CA at time t 
depends on the states of the connected neighbors at time t - 1. For a binary CA cell 
with a neighborhood size of N, there are 22 N possible functions. Table 1 illustrates the 

numbers involved. As Table 1 shows, the universe of possible functions increases 
extremely rapidly as the number of neighbors N grows. 
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Neighborhood size N 


# of Possible Functions 


1 


4 


2 


16 


3 


256 


4 


65,536 


5 


4,294,967,296 


6 


1.84 x 10 iy 


7 


3.4 x 10 JS 



Table 1 



It is theoretically possible to exhaustively search for viable CA-based RNG. 
However, in reality, the exhaustive search may be conducted for a small 
neighborhood size due to the tremendous growth of the search space (truth tables). 
With modern state of the art computing, N=4 may be the practical limit for exhaustive 
searches. 

SUMMARY OF THE INVENTION 

In a first aspect of the present invention, an embodiment of a method to reduce 
a search space for determining viable cellular automata based random number 
generators (CA-based RNGs) may include counting number of Is and Os of outputs of 
a truth table for a candidate CA-based RNG and counting number of Is and Os of 
inputs of the truth table for the candidate CA-based RNG. The method may also 
include accepting or rejecting the candidate CA-based RNG based on results of the 
counting steps. 

In a second aspect of the present invention, a system to reduce a search space 
for determining viable cellular automata based random number generator (CA-based 
RNGs) may include a truth-table-counting-module counting number of Is and Os of 
outputs of a truth table for a candidate CA-based RNG. The truth-table-counting 
module may also count number of Is and Os of inputs of the truth table for the 
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candidate CA-based RNG. The system may also include a prescreening-module 
accepting or rejecting the candidate CA-based RNG based on an output or outputs of 
the truth-table-counting module. 

In a third aspect of the present invention, computer readable medium may 
have embedded a software comprising a set of instructions for performing a method to 
reduce a search space for determining viable cellular automata based random number 
generator (CA-based RNGs). The embedded method may include counting number 
of Is and Os of outputs of a truth table for a candidate CA-based RNG and counting 
number of Is and Os of inputs of the truth table for the candidate CA-based RNG. 
The method may also include accepting or rejecting the candidate CA-based RNG 
based on results of the counting steps. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Features of the present invention will become apparent to those skilled in the 
art from the following description with reference to the drawings, in which: 

Figure 1A illustrates an exemplary truth table for a four-input cellular 
automata cell and the naming notation for the cellular automata; 

Figure IB illustrates an exemplary implementation of a cell of a cellular 
automata; 

Figure 1C illustrates an exemplary notation, a relative displacement notation, 
which provides a connection information of a CA cell; 

Figure ID illustrates an exemplary cellular automata showing the relationship 
between the relative displacement notation and the interconnection topology; 

Figure 2 illustrates an exemplary method to prescreen a candidate CA-based 
RNG; and 
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Figure 3 illustrates a block diagram of an exemplary system to prescreen a 
candidate CA-based RNG. 

DETAILED DESCRIPTION 

For simplicity and illustrative purposes, the principles of the present invention 
are described by referring mainly to exemplary embodiments thereof. However, one 
of ordinary skill in the art would readily recognize that the same principles are equally 
applicable to many situations in which random numbers generators are determined. 

High quality random numbers generators (RNGs) that pass stringent statistical 
tests may be implemented with cellular automata (CA). The basis of each cell is a 
logic function, which can be described by a truth table such as shown in Figure 1 A. It 
is also discussed above that the number of binary logic truth tables with N-inputs is 

22 N . As shown in Table 1, for N=4, the number of truth tables is 65,536. When N=5, 

the number of truth tables for a particular topology grows to over 4 billion. 

To put this into perspective, assume that viable CA-based RNGs with N=5 are 
being searched. The simplest instance is where the CA-based RNG has identical- 
function cells, i.e., the truth table is identical for all cells for the CA. In this instance, 
for a given topology, there are over 4 billion candidate RNGs, and each candidate 
RNG must be tested and evaluated. Depending of the length of the random number 
desired, the testing time will correspondingly increase. For example, desired length 
of the random may be 32 bits, 64 bits, etc. This process must be repeated for all 
possible topologies. As the numbers show, when searching for new random number 
generator implementations, reducing the search space is greatly desirable. 

After performing exhaustive searches on neighborhood size of 4 CA-based 
RNGs, the inventors of the present invention have discovered that the CA-based 
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RNGs that pass the battery of stringent random number tests (such as the DIEHARD 
suite of tests) all have common characteristics regarding their functions as represented 
by their truth tables. 

First, the number of Is in the output column was typically equal to the number 
5 of Os, i.e., each count was 8. Second, the number of Is and Os in the input 
contributing to output a 0 were typically equal as well. Similarly, the number of Is 
and Os in the input contributing to output a 1 were typically equal. This is clarified by 
the example below. 



Assume that a truth table is as follows (CA21530): 
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10 

For the CA21530, there are eight Is and eight Os in the output. Also, for all 
combination of inputs (dg, <U, d 2 , di) contributing to output a 0, there are sixteen Is 
and sixteen Os in the inputs. In a similar manner, it is also seen that there are sixteen 
Is and sixteen Os in the input contributing to output a 1. This indicates that the 

15 CA21530 is good candidate to pass the battery of random number tests, and thus 
passes the prescreening process. However, majority of the truth tables do not exhibit 
these characteristics and thus would not pass the prescreening process. This reduces 
the search space considerably. 

Figure 2 illustrates an exemplary method 200 to reduce a search space to 

20 determine viable CA-based RNGs. More specifically the method 200 prescreens a 
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candidate CA-based RNG. In step 210, Is and 0s of the truth table output of the 
candidate CA-based RNG may be counted. In step 220, if the difference in the output 
counts is less than or equal to a predetermined output difference threshold, then the 
method 200 may proceed to step 230. Otherwise, the method 200 may proceed to 
5 step 280 indicating that the particular candidate RNG has failed the prescreen process. 

The predetermined output difference threshold may be zero indicating that 
there must be equal number of Is and 0s. However, it is within the scope of the 
invention that strict adherence to equal number of Is and 0s may not be necessary, 
especially as the neighborhood size N grows larger. Thus, if N is 5 or greater, then 
' f J 10 perhaps a count difference of 2 or even 4 may be tolerated. Note this predetermined 

«! j output difference threshold is a parameter that may be set. 

-•4 

i,p In step 230, the method 200 counts the Is and 0s of the inputs in the truth table 

m 

i,H that generate Is as outputs. In step 240, if the difference in the input count is less than 

^ or equal to a predetermined Is input difference threshold, then the method 200 

15 proceeds to step 250. Else, the method 200 proceeds to step 280. Again, the 
predetermined Is difference threshold may be set to be greater than 0. 

In step 250, the method 200 counts the Is and 0s of the inputs in the truth table 
that generate 0s as outputs. In step 260, if the difference in the input count is less than 
or equal to a predetermined 0s input difference threshold, then the method 200 
20 proceeds to step 270 indicating that the candidate RNG has passed the prescreening 
process. Else, the method 200 proceeds to step 280. As before, the predetermined 0s 
difference threshold may be set to be greater than 0. 

Note that the steps of the method 200 may be modified and achieve a similar 
result. The steps may be modified, deleted or other steps may be added and still be 
25 within the scope of the invention. 



Q 
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The following example demonstrates how the screening process described 
above may reduce the search space. For a neighborhood size of 5 (each truth table 
has 32 entries), exhaustive search would require over 4 billion candidate RNGs to be 
evaluated for each given topology. However, if a strict equality of output counts is 
enforced, the number of candidate RNGs having sixteen Is and sixteen Os in the 
output is reduced to 601,080,390. In addition, if a strict equality of input counts is 
enforced, then the number of candidate RNGs is further reduced to 36,497,130. Thus 
from the original search space of 4,294,967,296, the search space is reduced by a 
factor of over 100 - a reduction of over two orders of magnitude. 

The method 200 may exist in a variety of forms both active and inactive. For 
example, they may exist as software program(s) comprised of program instructions in 
source code, object code, executable code or other formats. Any of the above may be 
embodied on a computer readable medium, which include storage devices and signals, 
in compressed or uncompressed form. Exemplary computer readable storage devices 
include conventional computer system RAM (random access memory), ROM (read 
only memory), EPROM (erasable, programmable ROM), EEPROM (electrically 
erasable, programmable ROM), flash memory, and magnetic or optical disks or tapes. 
Exemplary computer readable signals, whether modulated using a carrier or not, are 
signals that a computer system hosting or running the computer program may be 
configured to access, including signals downloaded through the Internet or other 
networks. Concrete examples of the foregoing include distribution of the program(s) 
on a CD ROM or via Internet download. In a sense, the Internet itself, as an abstract 
entity, is a computer readable medium. The same is true of computer networks in 
general. 
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Figure 4 illustrates a block diagram of an exemplary system 400 to prescreen a 
candidate CA-based RNG. As shown, the system may include a truth-table-counting- 
module 410 counting the outputs and the inputs of the truth table of the candidate CA- 
based RNG. The output counting may be performed by an output-counting-module 
412 and the input counting may be performed by an input-counting-module 414. The 
system 400 may also include a prescreening-module 420 which accepts or rejects the 
candidate CA-based RNG based on the results outputted by the truth-table-counting- 
module 410. 

While the invention has been described with reference to the exemplary 
embodiments thereof, those skilled in the art will be able to make various 
modifications to the described embodiments of the invention without departing from 
the true spirit and scope of the invention. 

The terms and descriptions used herein are set forth by way of illustration only 
and are not meant as limitations. In particular, although the method of the present 
invention has been described by examples, the steps of the method may be performed 
in a different order than illustrated or simultaneously. Those skilled in the art will 
recognize that these and other variations are possible within the spirit and scope of 
the invention as defined in the following claims and their equivalents. 
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